{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# p27: Solve KdV equation"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Solve the KdV equation using FFT\n",
    "$$\n",
    "u_t + u u_x + u_{xxx} = 0, \\qquad x \\in [-\\pi,\\pi]\n",
    "$$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "%config InlineBackend.figure_format='svg'\n",
    "from mpl_toolkits.mplot3d import Axes3D\n",
    "from matplotlib.collections import LineCollection\n",
    "from numpy import pi,cosh,exp,round,zeros,arange,real\n",
    "from numpy.fft import fft,ifft\n",
    "from matplotlib.pyplot import figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?>\n",
       "<!DOCTYPE svg PUBLIC \"-//W3C//DTD SVG 1.1//EN\"\n",
       "  \"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd\">\n",
       "<!-- Created with matplotlib (http://matplotlib.org/) -->\n",
       "<svg height=\"349pt\" version=\"1.1\" viewBox=\"0 0 683 349\" width=\"683pt\" xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\">\n",
       " <defs>\n",
       "  <style type=\"text/css\">\n",
       "*{stroke-linecap:butt;stroke-linejoin:round;stroke-miterlimit:100000;}\n",
       "  </style>\n",
       " </defs>\n",
       " <g id=\"figure_1\">\n",
       "  <g id=\"patch_1\">\n",
       "   <path d=\"M 0 349.2 \n",
       "L 684 349.2 \n",
       "L 684 0 \n",
       "L 0 0 \n",
       "z\n",
       "\" style=\"fill:none;\"/>\n",
       "  </g>\n",
       "  <g id=\"patch_2\">\n",
       "   <path d=\"M 7.2 342 \n",
       "L 676.8 342 \n",
       "L 676.8 7.2 \n",
       "L 7.2 7.2 \n",
       "z\n",
       "\" style=\"fill:#ffffff;\"/>\n",
       "  </g>\n",
       "  <g id=\"pane3d_1\">\n",
       "   <g id=\"patch_3\">\n",
       "    <path d=\"M 601.315913 207.406815 \n",
       "L 394.562923 63.944023 \n",
       "L 399.206484 18.565608 \n",
       "L 628.440907 176.867094 \n",
       "\" style=\"fill:#f2f2f2;opacity:0.5;stroke:#f2f2f2;stroke-linejoin:miter;\"/>\n",
       "   </g>\n",
       "  </g>\n",
       "  <g id=\"pane3d_2\">\n",
       "   <g id=\"patch_4\">\n",
       "    <path d=\"M 101.836618 163.974196 \n",
       "L 394.562923 63.944023 \n",
       "L 399.206484 18.565608 \n",
       "L 74.952168 128.891719 \n",
       "\" style=\"fill:#e6e6e6;opacity:0.5;stroke:#e6e6e6;stroke-linejoin:miter;\"/>\n",
       "   </g>\n",
       "  </g>\n",
       "  <g id=\"pane3d_3\">\n",
       "   <g id=\"patch_5\">\n",
       "    <path d=\"M 306.478995 310.296912 \n",
       "L 601.315913 207.406815 \n",
       "L 394.562923 63.944023 \n",
       "L 101.836618 163.974196 \n",
       "\" style=\"fill:#ececec;opacity:0.5;stroke:#ececec;stroke-linejoin:miter;\"/>\n",
       "   </g>\n",
       "  </g>\n",
       "  <g id=\"axis3d_1\">\n",
       "   <g id=\"line2d_1\">\n",
       "    <path d=\"M 601.315913 207.406815 \n",
       "L 306.478995 310.296912 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-width:0.75;\"/>\n",
       "   </g>\n",
       "   <g id=\"text_1\">\n",
       "    <!-- x -->\n",
       "    <defs>\n",
       "     <path d=\"M 54.890625 54.6875 \n",
       "L 35.109375 28.078125 \n",
       "L 55.90625 0 \n",
       "L 45.3125 0 \n",
       "L 29.390625 21.484375 \n",
       "L 13.484375 0 \n",
       "L 2.875 0 \n",
       "L 24.125 28.609375 \n",
       "L 4.6875 54.6875 \n",
       "L 15.28125 54.6875 \n",
       "L 29.78125 35.203125 \n",
       "L 44.28125 54.6875 \n",
       "z\n",
       "\" id=\"BitstreamVeraSans-Roman-78\"/>\n",
       "    </defs>\n",
       "    <g transform=\"translate(471.234008695 278.074392759)scale(0.1 -0.1)\">\n",
       "     <use xlink:href=\"#BitstreamVeraSans-Roman-78\"/>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"Line3DCollection_1\">\n",
       "    <path d=\"M 318.871623 305.97222 \n",
       "L 114.138893 159.770273 \n",
       "L 88.593189 124.250421 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 364.262285 290.132109 \n",
       "L 159.199846 144.372084 \n",
       "L 138.547425 107.253711 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 409.514882 274.340178 \n",
       "L 204.125658 129.020076 \n",
       "L 188.335742 90.313455 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 454.630042 258.59621 \n",
       "L 248.916935 113.714041 \n",
       "L 237.958967 73.42937 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 499.608391 242.899985 \n",
       "L 293.574281 98.453772 \n",
       "L 287.417918 56.601178 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 544.450549 227.251286 \n",
       "L 338.098295 83.239066 \n",
       "L 336.713411 39.828603 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 589.157134 211.649899 \n",
       "L 382.489574 68.069717 \n",
       "L 385.846254 23.111368 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "   </g>\n",
       "   <g id=\"xtick_1\">\n",
       "    <g id=\"line2d_2\">\n",
       "     <path d=\"M 317.210596 304.786062 \n",
       "L 322.194814 308.345348 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_2\">\n",
       "     <!-- −3 -->\n",
       "     <defs>\n",
       "      <path d=\"M 10.59375 35.5 \n",
       "L 73.1875 35.5 \n",
       "L 73.1875 27.203125 \n",
       "L 10.59375 27.203125 \n",
       "z\n",
       "\" id=\"BitstreamVeraSans-Roman-2212\"/>\n",
       "      <path d=\"M 40.578125 39.3125 \n",
       "Q 47.65625 37.796875 51.625 33 \n",
       "Q 55.609375 28.21875 55.609375 21.1875 \n",
       "Q 55.609375 10.40625 48.1875 4.484375 \n",
       "Q 40.765625 -1.421875 27.09375 -1.421875 \n",
       "Q 22.515625 -1.421875 17.65625 -0.515625 \n",
       "Q 12.796875 0.390625 7.625 2.203125 \n",
       "L 7.625 11.71875 \n",
       "Q 11.71875 9.328125 16.59375 8.109375 \n",
       "Q 21.484375 6.890625 26.8125 6.890625 \n",
       "Q 36.078125 6.890625 40.9375 10.546875 \n",
       "Q 45.796875 14.203125 45.796875 21.1875 \n",
       "Q 45.796875 27.640625 41.28125 31.265625 \n",
       "Q 36.765625 34.90625 28.71875 34.90625 \n",
       "L 20.21875 34.90625 \n",
       "L 20.21875 43.015625 \n",
       "L 29.109375 43.015625 \n",
       "Q 36.375 43.015625 40.234375 45.921875 \n",
       "Q 44.09375 48.828125 44.09375 54.296875 \n",
       "Q 44.09375 59.90625 40.109375 62.90625 \n",
       "Q 36.140625 65.921875 28.71875 65.921875 \n",
       "Q 24.65625 65.921875 20.015625 65.03125 \n",
       "Q 15.375 64.15625 9.8125 62.3125 \n",
       "L 9.8125 71.09375 \n",
       "Q 15.4375 72.65625 20.34375 73.4375 \n",
       "Q 25.25 74.21875 29.59375 74.21875 \n",
       "Q 40.828125 74.21875 47.359375 69.109375 \n",
       "Q 53.90625 64.015625 53.90625 55.328125 \n",
       "Q 53.90625 49.265625 50.4375 45.09375 \n",
       "Q 46.96875 40.921875 40.578125 39.3125 \n",
       "\" id=\"BitstreamVeraSans-Roman-33\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(321.135617128 321.174086091)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-2212\"/>\n",
       "      <use x=\"83.7890625\" xlink:href=\"#BitstreamVeraSans-Roman-33\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_2\">\n",
       "    <g id=\"line2d_3\">\n",
       "     <path d=\"M 362.598619 288.949561 \n",
       "L 367.590756 292.498012 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_3\">\n",
       "     <!-- −2 -->\n",
       "     <defs>\n",
       "      <path d=\"M 19.1875 8.296875 \n",
       "L 53.609375 8.296875 \n",
       "L 53.609375 0 \n",
       "L 7.328125 0 \n",
       "L 7.328125 8.296875 \n",
       "Q 12.9375 14.109375 22.625 23.890625 \n",
       "Q 32.328125 33.6875 34.8125 36.53125 \n",
       "Q 39.546875 41.84375 41.421875 45.53125 \n",
       "Q 43.3125 49.21875 43.3125 52.78125 \n",
       "Q 43.3125 58.59375 39.234375 62.25 \n",
       "Q 35.15625 65.921875 28.609375 65.921875 \n",
       "Q 23.96875 65.921875 18.8125 64.3125 \n",
       "Q 13.671875 62.703125 7.8125 59.421875 \n",
       "L 7.8125 69.390625 \n",
       "Q 13.765625 71.78125 18.9375 73 \n",
       "Q 24.125 74.21875 28.421875 74.21875 \n",
       "Q 39.75 74.21875 46.484375 68.546875 \n",
       "Q 53.21875 62.890625 53.21875 53.421875 \n",
       "Q 53.21875 48.921875 51.53125 44.890625 \n",
       "Q 49.859375 40.875 45.40625 35.40625 \n",
       "Q 44.1875 33.984375 37.640625 27.21875 \n",
       "Q 31.109375 20.453125 19.1875 8.296875 \n",
       "\" id=\"BitstreamVeraSans-Roman-32\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(366.336717769 305.383487982)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-2212\"/>\n",
       "      <use x=\"83.7890625\" xlink:href=\"#BitstreamVeraSans-Roman-32\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_3\">\n",
       "    <g id=\"line2d_4\">\n",
       "     <path d=\"M 407.8486 273.161225 \n",
       "L 412.848584 276.69889 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_4\">\n",
       "     <!-- −1 -->\n",
       "     <defs>\n",
       "      <path d=\"M 12.40625 8.296875 \n",
       "L 28.515625 8.296875 \n",
       "L 28.515625 63.921875 \n",
       "L 10.984375 60.40625 \n",
       "L 10.984375 69.390625 \n",
       "L 28.421875 72.90625 \n",
       "L 38.28125 72.90625 \n",
       "L 38.28125 8.296875 \n",
       "L 54.390625 8.296875 \n",
       "L 54.390625 0 \n",
       "L 12.40625 0 \n",
       "z\n",
       "\" id=\"BitstreamVeraSans-Roman-31\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(411.400858297 289.640735656)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-2212\"/>\n",
       "      <use x=\"83.7890625\" xlink:href=\"#BitstreamVeraSans-Roman-31\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_4\">\n",
       "    <g id=\"line2d_5\">\n",
       "     <path d=\"M 452.961167 257.420834 \n",
       "L 457.968929 260.947762 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_5\">\n",
       "     <!-- 0 -->\n",
       "     <defs>\n",
       "      <path d=\"M 31.78125 66.40625 \n",
       "Q 24.171875 66.40625 20.328125 58.90625 \n",
       "Q 16.5 51.421875 16.5 36.375 \n",
       "Q 16.5 21.390625 20.328125 13.890625 \n",
       "Q 24.171875 6.390625 31.78125 6.390625 \n",
       "Q 39.453125 6.390625 43.28125 13.890625 \n",
       "Q 47.125 21.390625 47.125 36.375 \n",
       "Q 47.125 51.421875 43.28125 58.90625 \n",
       "Q 39.453125 66.40625 31.78125 66.40625 \n",
       "M 31.78125 74.21875 \n",
       "Q 44.046875 74.21875 50.515625 64.515625 \n",
       "Q 56.984375 54.828125 56.984375 36.375 \n",
       "Q 56.984375 17.96875 50.515625 8.265625 \n",
       "Q 44.046875 -1.421875 31.78125 -1.421875 \n",
       "Q 19.53125 -1.421875 13.0625 8.265625 \n",
       "Q 6.59375 17.96875 6.59375 36.375 \n",
       "Q 6.59375 54.828125 13.0625 64.515625 \n",
       "Q 19.53125 74.21875 31.78125 74.21875 \n",
       "\" id=\"BitstreamVeraSans-Roman-30\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(460.518504009 273.945611982)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_5\">\n",
       "    <g id=\"line2d_6\">\n",
       "     <path d=\"M 497.936947 241.72817 \n",
       "L 502.952416 245.244411 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_6\">\n",
       "     <!-- 1 -->\n",
       "     <g transform=\"translate(505.310585194 258.297901142)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-31\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_6\">\n",
       "    <g id=\"line2d_7\">\n",
       "     <path d=\"M 542.776559 226.083017 \n",
       "L 547.799665 229.588618 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_7\">\n",
       "     <!-- 2 -->\n",
       "     <g transform=\"translate(549.967559667 242.697388616)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-32\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_7\">\n",
       "    <g id=\"line2d_8\">\n",
       "     <path d=\"M 587.480622 210.485159 \n",
       "L 592.511296 213.980169 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_8\">\n",
       "     <!-- 3 -->\n",
       "     <g transform=\"translate(594.490037792 227.143861181)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-33\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "  </g>\n",
       "  <g id=\"axis3d_2\">\n",
       "   <g id=\"line2d_9\">\n",
       "    <path d=\"M 101.836618 163.974196 \n",
       "L 306.478995 310.296912 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-width:0.75;\"/>\n",
       "   </g>\n",
       "   <g id=\"text_9\">\n",
       "    <!-- t -->\n",
       "    <defs>\n",
       "     <path d=\"M 18.3125 70.21875 \n",
       "L 18.3125 54.6875 \n",
       "L 36.8125 54.6875 \n",
       "L 36.8125 47.703125 \n",
       "L 18.3125 47.703125 \n",
       "L 18.3125 18.015625 \n",
       "Q 18.3125 11.328125 20.140625 9.421875 \n",
       "Q 21.96875 7.515625 27.59375 7.515625 \n",
       "L 36.8125 7.515625 \n",
       "L 36.8125 0 \n",
       "L 27.59375 0 \n",
       "Q 17.1875 0 13.234375 3.875 \n",
       "Q 9.28125 7.765625 9.28125 18.015625 \n",
       "L 9.28125 47.703125 \n",
       "L 2.6875 47.703125 \n",
       "L 2.6875 54.6875 \n",
       "L 9.28125 54.6875 \n",
       "L 9.28125 70.21875 \n",
       "z\n",
       "\" id=\"BitstreamVeraSans-Roman-74\"/>\n",
       "    </defs>\n",
       "    <g transform=\"translate(173.663324104 251.912777812)scale(0.1 -0.1)\">\n",
       "     <use xlink:href=\"#BitstreamVeraSans-Roman-74\"/>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"Line3DCollection_2\">\n",
       "    <path d=\"M 623.785904 173.652507 \n",
       "L 597.123633 204.497855 \n",
       "L 302.328946 307.329558 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 586.650518 148.008075 \n",
       "L 563.670448 181.285193 \n",
       "L 269.21359 283.651525 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 549.700179 122.491428 \n",
       "L 530.367738 158.176943 \n",
       "L 236.24867 260.081057 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 512.933506 97.101617 \n",
       "L 497.214491 135.172403 \n",
       "L 203.433163 236.617422 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 476.349135 71.837696 \n",
       "L 464.209702 112.270876 \n",
       "L 170.766057 213.259896 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 439.945712 46.698731 \n",
       "L 431.352376 89.471672 \n",
       "L 138.246347 190.00776 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 403.721899 21.683799 \n",
       "L 398.641528 66.774106 \n",
       "L 105.873037 166.860303 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "   </g>\n",
       "   <g id=\"xtick_8\">\n",
       "    <g id=\"line2d_10\">\n",
       "     <path d=\"M 304.710663 306.498756 \n",
       "L 297.56437 308.99156 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_10\">\n",
       "     <!-- 0.000 -->\n",
       "     <defs>\n",
       "      <path d=\"M 10.6875 12.40625 \n",
       "L 21 12.40625 \n",
       "L 21 0 \n",
       "L 10.6875 0 \n",
       "z\n",
       "\" id=\"BitstreamVeraSans-Roman-2e\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(274.630384689 320.493905038)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-2e\"/>\n",
       "      <use x=\"95.41015625\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"159.033203125\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"222.65625\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_9\">\n",
       "    <g id=\"line2d_11\">\n",
       "     <path d=\"M 271.592524 282.824501 \n",
       "L 264.454583 285.305968 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_11\">\n",
       "     <!-- 0.001 -->\n",
       "     <g transform=\"translate(241.6793575 296.899051628)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-2e\"/>\n",
       "      <use x=\"95.41015625\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"159.033203125\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"222.65625\" xlink:href=\"#BitstreamVeraSans-Roman-31\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_10\">\n",
       "    <g id=\"line2d_12\">\n",
       "     <path d=\"M 238.624822 259.257786 \n",
       "L 231.49523 261.727992 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_12\">\n",
       "     <!-- 0.002 -->\n",
       "     <g transform=\"translate(208.877416019 273.410952263)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-2e\"/>\n",
       "      <use x=\"95.41015625\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"159.033203125\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"222.65625\" xlink:href=\"#BitstreamVeraSans-Roman-32\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_11\">\n",
       "    <g id=\"line2d_13\">\n",
       "     <path d=\"M 205.806535 235.797878 \n",
       "L 198.68529 238.2569 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_13\">\n",
       "     <!-- 0.003 -->\n",
       "     <g transform=\"translate(176.22355073 250.028884071)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-2e\"/>\n",
       "      <use x=\"95.41015625\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"159.033203125\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"222.65625\" xlink:href=\"#BitstreamVeraSans-Roman-33\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_12\">\n",
       "    <g id=\"line2d_14\">\n",
       "     <path d=\"M 173.136649 212.444054 \n",
       "L 166.023748 214.891968 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_14\">\n",
       "     <!-- 0.004 -->\n",
       "     <defs>\n",
       "      <path d=\"M 37.796875 64.3125 \n",
       "L 12.890625 25.390625 \n",
       "L 37.796875 25.390625 \n",
       "z\n",
       "M 35.203125 72.90625 \n",
       "L 47.609375 72.90625 \n",
       "L 47.609375 25.390625 \n",
       "L 58.015625 25.390625 \n",
       "L 58.015625 17.1875 \n",
       "L 47.609375 17.1875 \n",
       "L 47.609375 0 \n",
       "L 37.796875 0 \n",
       "L 37.796875 17.1875 \n",
       "L 4.890625 17.1875 \n",
       "L 4.890625 26.703125 \n",
       "z\n",
       "\" id=\"BitstreamVeraSans-Roman-34\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(143.716761211 226.75213069)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-2e\"/>\n",
       "      <use x=\"95.41015625\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"159.033203125\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"222.65625\" xlink:href=\"#BitstreamVeraSans-Roman-34\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_13\">\n",
       "    <g id=\"line2d_15\">\n",
       "     <path d=\"M 140.614159 189.195595 \n",
       "L 133.5096 191.632476 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_15\">\n",
       "     <!-- 0.005 -->\n",
       "     <defs>\n",
       "      <path d=\"M 10.796875 72.90625 \n",
       "L 49.515625 72.90625 \n",
       "L 49.515625 64.59375 \n",
       "L 19.828125 64.59375 \n",
       "L 19.828125 46.734375 \n",
       "Q 21.96875 47.46875 24.109375 47.828125 \n",
       "Q 26.265625 48.1875 28.421875 48.1875 \n",
       "Q 40.625 48.1875 47.75 41.5 \n",
       "Q 54.890625 34.8125 54.890625 23.390625 \n",
       "Q 54.890625 11.625 47.5625 5.09375 \n",
       "Q 40.234375 -1.421875 26.90625 -1.421875 \n",
       "Q 22.3125 -1.421875 17.546875 -0.640625 \n",
       "Q 12.796875 0.140625 7.71875 1.703125 \n",
       "L 7.71875 11.625 \n",
       "Q 12.109375 9.234375 16.796875 8.0625 \n",
       "Q 21.484375 6.890625 26.703125 6.890625 \n",
       "Q 35.15625 6.890625 40.078125 11.328125 \n",
       "Q 45.015625 15.765625 45.015625 23.390625 \n",
       "Q 45.015625 31 40.078125 35.4375 \n",
       "Q 35.15625 39.890625 26.703125 39.890625 \n",
       "Q 22.75 39.890625 18.8125 39.015625 \n",
       "Q 14.890625 38.140625 10.796875 36.28125 \n",
       "z\n",
       "\" id=\"BitstreamVeraSans-Roman-35\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(111.356056033 203.579982201)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-2e\"/>\n",
       "      <use x=\"95.41015625\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"159.033203125\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"222.65625\" xlink:href=\"#BitstreamVeraSans-Roman-35\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_14\">\n",
       "    <g id=\"line2d_16\">\n",
       "     <path d=\"M 108.23807 166.05179 \n",
       "L 101.141851 168.477712 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_16\">\n",
       "     <!-- 0.006 -->\n",
       "     <defs>\n",
       "      <path d=\"M 33.015625 40.375 \n",
       "Q 26.375 40.375 22.484375 35.828125 \n",
       "Q 18.609375 31.296875 18.609375 23.390625 \n",
       "Q 18.609375 15.53125 22.484375 10.953125 \n",
       "Q 26.375 6.390625 33.015625 6.390625 \n",
       "Q 39.65625 6.390625 43.53125 10.953125 \n",
       "Q 47.40625 15.53125 47.40625 23.390625 \n",
       "Q 47.40625 31.296875 43.53125 35.828125 \n",
       "Q 39.65625 40.375 33.015625 40.375 \n",
       "M 52.59375 71.296875 \n",
       "L 52.59375 62.3125 \n",
       "Q 48.875 64.0625 45.09375 64.984375 \n",
       "Q 41.3125 65.921875 37.59375 65.921875 \n",
       "Q 27.828125 65.921875 22.671875 59.328125 \n",
       "Q 17.53125 52.734375 16.796875 39.40625 \n",
       "Q 19.671875 43.65625 24.015625 45.921875 \n",
       "Q 28.375 48.1875 33.59375 48.1875 \n",
       "Q 44.578125 48.1875 50.953125 41.515625 \n",
       "Q 57.328125 34.859375 57.328125 23.390625 \n",
       "Q 57.328125 12.15625 50.6875 5.359375 \n",
       "Q 44.046875 -1.421875 33.015625 -1.421875 \n",
       "Q 20.359375 -1.421875 13.671875 8.265625 \n",
       "Q 6.984375 17.96875 6.984375 36.375 \n",
       "Q 6.984375 53.65625 15.1875 63.9375 \n",
       "Q 23.390625 74.21875 37.203125 74.21875 \n",
       "Q 40.921875 74.21875 44.703125 73.484375 \n",
       "Q 48.484375 72.75 52.59375 71.296875 \n",
       "\" id=\"BitstreamVeraSans-Roman-36\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(79.1404526564 180.511735047)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-2e\"/>\n",
       "      <use x=\"95.41015625\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"159.033203125\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"222.65625\" xlink:href=\"#BitstreamVeraSans-Roman-36\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "  </g>\n",
       "  <g id=\"axis3d_3\">\n",
       "   <g id=\"line2d_17\">\n",
       "    <path d=\"M 101.836618 163.974196 \n",
       "L 74.952168 128.891719 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-width:0.75;\"/>\n",
       "   </g>\n",
       "   <g id=\"text_17\">\n",
       "    <!-- u -->\n",
       "    <defs>\n",
       "     <path d=\"M 8.5 21.578125 \n",
       "L 8.5 54.6875 \n",
       "L 17.484375 54.6875 \n",
       "L 17.484375 21.921875 \n",
       "Q 17.484375 14.15625 20.5 10.265625 \n",
       "Q 23.53125 6.390625 29.59375 6.390625 \n",
       "Q 36.859375 6.390625 41.078125 11.03125 \n",
       "Q 45.3125 15.671875 45.3125 23.6875 \n",
       "L 45.3125 54.6875 \n",
       "L 54.296875 54.6875 \n",
       "L 54.296875 0 \n",
       "L 45.3125 0 \n",
       "L 45.3125 8.40625 \n",
       "Q 42.046875 3.421875 37.71875 1 \n",
       "Q 33.40625 -1.421875 27.6875 -1.421875 \n",
       "Q 18.265625 -1.421875 13.375 4.4375 \n",
       "Q 8.5 10.296875 8.5 21.578125 \n",
       "\" id=\"BitstreamVeraSans-Roman-75\"/>\n",
       "    </defs>\n",
       "    <g transform=\"translate(33.9716154336 145.580250008)scale(0.1 -0.1)\">\n",
       "     <use xlink:href=\"#BitstreamVeraSans-Roman-75\"/>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"Line3DCollection_3\">\n",
       "    <path d=\"M 101.350338 163.339633 \n",
       "L 394.647001 63.122382 \n",
       "L 601.806323 206.854668 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 95.363414 155.527082 \n",
       "L 395.681886 53.009141 \n",
       "L 607.844769 200.05605 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 89.082344 147.33069 \n",
       "L 396.767096 42.404115 \n",
       "L 614.181213 192.921921 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 82.484906 138.721458 \n",
       "L 397.906392 31.270539 \n",
       "L 620.838271 185.426815 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "    <path d=\"M 75.546579 129.667387 \n",
       "L 399.103921 19.567887 \n",
       "L 627.84091 177.542625 \n",
       "\" style=\"fill:none;stroke:#e6e6e6;\"/>\n",
       "   </g>\n",
       "   <g id=\"xtick_15\">\n",
       "    <g id=\"line2d_18\">\n",
       "     <path d=\"M 103.719677 162.530047 \n",
       "L 96.610539 164.959186 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_18\">\n",
       "     <!-- 0 -->\n",
       "     <g transform=\"translate(75.3105287281 168.950500588)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_16\">\n",
       "    <g id=\"line2d_19\">\n",
       "     <path d=\"M 97.79004 154.698718 \n",
       "L 90.508983 157.18421 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_19\">\n",
       "     <!-- 500 -->\n",
       "     <g transform=\"translate(62.4143297566 161.091854038)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-35\"/>\n",
       "      <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"127.24609375\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_17\">\n",
       "    <g id=\"line2d_20\">\n",
       "     <path d=\"M 91.569098 146.482658 \n",
       "L 84.107601 149.027175 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_20\">\n",
       "     <!-- 1000 -->\n",
       "     <g transform=\"translate(52.3784371462 152.847180731)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-31\"/>\n",
       "      <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"127.24609375\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"190.869140625\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_18\">\n",
       "    <g id=\"line2d_21\">\n",
       "     <path d=\"M 85.034841 137.852801 \n",
       "L 77.383736 140.459213 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_21\">\n",
       "     <!-- 1500 -->\n",
       "     <g transform=\"translate(45.178607727 144.187321304)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-31\"/>\n",
       "      <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-35\"/>\n",
       "      <use x=\"127.24609375\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"190.869140625\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"xtick_19\">\n",
       "    <g id=\"line2d_22\">\n",
       "     <path d=\"M 78.16299 128.77708 \n",
       "L 70.312389 131.448468 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "    </g>\n",
       "    <g id=\"text_22\">\n",
       "     <!-- 2000 -->\n",
       "     <g transform=\"translate(37.6068435818 135.080103717)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-32\"/>\n",
       "      <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"127.24609375\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      <use x=\"190.869140625\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "  </g>\n",
       "  <g id=\"axes_1\">\n",
       "   <g id=\"Line3DCollection_4\">\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 308.197926 304.888958 \n",
       "L 348.40038 290.754243 \n",
       "L 350.631445 289.709214 \n",
       "L 351.74973 289.003821 \n",
       "L 352.873779 288.041671 \n",
       "L 354.010621 286.633044 \n",
       "L 355.173427 284.489965 \n",
       "L 356.383172 281.265126 \n",
       "L 357.662763 276.79084 \n",
       "L 359.010567 271.666149 \n",
       "L 360.357489 267.780681 \n",
       "L 361.580007 267.305952 \n",
       "L 362.619753 270.296939 \n",
       "L 364.456354 278.363712 \n",
       "L 365.408889 280.864611 \n",
       "L 366.411468 282.274749 \n",
       "L 367.454308 282.9372 \n",
       "L 368.524738 283.142148 \n",
       "L 369.612645 283.082628 \n",
       "L 371.815943 282.584767 \n",
       "L 376.258879 281.137787 \n",
       "L 388.511066 276.802651 \n",
       "L 391.8814 275.467315 \n",
       "L 394.165642 274.402865 \n",
       "L 396.520658 273.026837 \n",
       "L 397.744112 272.144687 \n",
       "L 399.011616 271.083215 \n",
       "L 401.708482 268.336812 \n",
       "L 404.558002 265.104512 \n",
       "L 405.939463 263.736093 \n",
       "L 407.208032 262.858714 \n",
       "L 408.320781 262.61493 \n",
       "L 409.280345 262.960244 \n",
       "L 410.132152 263.690217 \n",
       "L 411.752507 265.341615 \n",
       "L 412.605066 265.952457 \n",
       "L 413.507621 266.352076 \n",
       "L 414.458821 266.555463 \n",
       "L 415.451419 266.596978 \n",
       "L 417.526403 266.339347 \n",
       "L 419.674416 265.812856 \n",
       "L 424.062265 264.433869 \n",
       "L 437.357016 259.836404 \n",
       "L 590.65844 206.363448 \n",
       "L 590.65844 206.363448 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 300.283675 299.231792 \n",
       "L 339.362284 285.617886 \n",
       "L 347.16313 282.776703 \n",
       "L 349.391052 281.735757 \n",
       "L 350.506538 281.030748 \n",
       "L 351.62602 280.075875 \n",
       "L 352.755277 278.650395 \n",
       "L 353.905093 276.508515 \n",
       "L 355.094356 273.24733 \n",
       "L 356.345048 268.686116 \n",
       "L 357.66089 263.393063 \n",
       "L 358.98955 259.210058 \n",
       "L 360.219755 258.49109 \n",
       "L 361.2879 261.442772 \n",
       "L 363.177332 269.90805 \n",
       "L 364.143711 272.592182 \n",
       "L 365.153237 274.134619 \n",
       "L 366.199379 274.877774 \n",
       "L 367.271659 275.122409 \n",
       "L 368.359361 275.095636 \n",
       "L 370.562331 274.615709 \n",
       "L 373.890526 273.560846 \n",
       "L 383.912275 269.991879 \n",
       "L 387.283898 268.614044 \n",
       "L 389.571848 267.481437 \n",
       "L 391.935356 265.971316 \n",
       "L 393.167463 264.975002 \n",
       "L 394.440892 263.793047 \n",
       "L 397.140349 260.780956 \n",
       "L 399.926336 257.551902 \n",
       "L 401.237885 256.392301 \n",
       "L 402.424243 255.832465 \n",
       "L 403.466327 255.920718 \n",
       "L 404.383597 256.533683 \n",
       "L 406.924676 259.013048 \n",
       "L 407.824163 259.538029 \n",
       "L 408.766733 259.859039 \n",
       "L 409.753034 259.984325 \n",
       "L 411.814446 259.847342 \n",
       "L 413.954973 259.379178 \n",
       "L 418.337518 258.036154 \n",
       "L 427.192809 254.997384 \n",
       "L 446.02188 248.442729 \n",
       "L 470.359636 239.956061 \n",
       "L 500.168063 229.572176 \n",
       "L 521.115901 222.264721 \n",
       "L 551.921868 211.531769 \n",
       "L 582.663444 200.820524 \n",
       "L 582.663444 200.820524 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 292.381218 293.591942 \n",
       "L 305.78353 288.909532 \n",
       "L 313.599935 286.196505 \n",
       "L 332.561074 279.594201 \n",
       "L 344.815463 275.269937 \n",
       "L 347.040546 274.344765 \n",
       "L 348.152797 273.783019 \n",
       "L 349.265515 273.06957 \n",
       "L 350.380785 272.119611 \n",
       "L 351.502745 270.693525 \n",
       "L 352.640536 268.536939 \n",
       "L 353.810152 265.257608 \n",
       "L 355.032907 260.612539 \n",
       "L 356.317135 255.148958 \n",
       "L 357.624964 250.682699 \n",
       "L 358.860167 249.692759 \n",
       "L 359.954494 252.594008 \n",
       "L 361.898003 261.442908 \n",
       "L 362.878653 264.322811 \n",
       "L 363.896252 265.994625 \n",
       "L 364.945816 266.827103 \n",
       "L 366.020199 267.111566 \n",
       "L 367.108111 267.116408 \n",
       "L 369.310921 266.654031 \n",
       "L 372.637503 265.604873 \n",
       "L 378.201862 263.616752 \n",
       "L 381.560279 262.278365 \n",
       "L 383.830397 261.185682 \n",
       "L 386.156022 259.766661 \n",
       "L 387.354463 258.858541 \n",
       "L 388.590026 257.750444 \n",
       "L 391.192792 254.879281 \n",
       "L 393.924721 251.500211 \n",
       "L 395.255909 250.092808 \n",
       "L 396.498718 249.199466 \n",
       "L 397.610681 249.003405 \n",
       "L 398.603194 249.400016 \n",
       "L 399.504534 250.210909 \n",
       "L 401.240444 251.995355 \n",
       "L 402.137862 252.671842 \n",
       "L 403.08025 253.101383 \n",
       "L 404.059552 253.334441 \n",
       "L 405.070521 253.405271 \n",
       "L 407.169691 253.170932 \n",
       "L 410.419033 252.346277 \n",
       "L 418.145845 249.777004 \n",
       "L 431.437451 245.154072 \n",
       "L 456.880664 236.298854 \n",
       "L 483.390223 227.059592 \n",
       "L 493.308801 223.619051 \n",
       "L 504.332524 219.777764 \n",
       "L 535.157219 209.035973 \n",
       "L 550.534891 203.689498 \n",
       "L 560.429538 200.23066 \n",
       "L 570.291335 196.8103 \n",
       "L 574.693162 195.262981 \n",
       "L 574.693162 195.262981 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 284.480796 287.936771 \n",
       "L 295.65131 284.053229 \n",
       "L 342.470603 267.724516 \n",
       "L 344.69471 266.874527 \n",
       "L 346.916625 265.842101 \n",
       "L 348.026672 265.124039 \n",
       "L 349.137964 264.173768 \n",
       "L 350.253115 262.753419 \n",
       "L 351.3795 260.581534 \n",
       "L 352.530583 257.286816 \n",
       "L 353.726251 252.570103 \n",
       "L 354.979121 246.935101 \n",
       "L 356.264084 242.198491 \n",
       "L 357.501307 240.916082 \n",
       "L 358.619798 243.75038 \n",
       "L 360.618294 252.967063 \n",
       "L 361.6138 256.055468 \n",
       "L 362.640338 257.855676 \n",
       "L 363.693984 258.778493 \n",
       "L 364.770069 259.113974 \n",
       "L 365.859319 259.136943 \n",
       "L 366.957132 258.978594 \n",
       "L 370.277918 258.013852 \n",
       "L 374.730284 256.390609 \n",
       "L 378.102024 254.900754 \n",
       "L 380.397196 253.564338 \n",
       "L 381.571681 252.716352 \n",
       "L 382.775717 251.67758 \n",
       "L 384.011601 250.449852 \n",
       "L 386.608576 247.293875 \n",
       "L 389.279608 243.931086 \n",
       "L 390.549925 242.761975 \n",
       "L 391.725842 242.206679 \n",
       "L 392.784783 242.350753 \n",
       "L 393.744467 243.039121 \n",
       "L 396.451434 245.720357 \n",
       "L 397.39374 246.2869 \n",
       "L 398.367203 246.651111 \n",
       "L 399.379539 246.787552 \n",
       "L 400.412339 246.801104 \n",
       "L 402.541501 246.479275 \n",
       "L 404.71324 245.928413 \n",
       "L 409.11102 244.547138 \n",
       "L 489.766918 216.516657 \n",
       "L 499.689247 213.052464 \n",
       "L 509.588664 209.620729 \n",
       "L 519.500038 206.158556 \n",
       "L 529.384756 202.732351 \n",
       "L 542.579317 198.127182 \n",
       "L 550.258471 195.461632 \n",
       "L 566.705503 189.745707 \n",
       "L 566.705503 189.745707 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 276.596091 282.30966 \n",
       "L 301.148814 273.765734 \n",
       "L 318.982139 267.568734 \n",
       "L 341.241356 259.79252 \n",
       "L 344.573289 258.459733 \n",
       "L 345.682367 257.906945 \n",
       "L 346.790082 257.196874 \n",
       "L 347.897473 256.238578 \n",
       "L 349.0063 254.824261 \n",
       "L 350.121891 252.646437 \n",
       "L 351.255471 249.328787 \n",
       "L 352.425016 244.556665 \n",
       "L 353.646667 238.75738 \n",
       "L 354.907344 233.74832 \n",
       "L 356.143127 232.177028 \n",
       "L 357.284181 234.905593 \n",
       "L 359.338063 244.480927 \n",
       "L 360.349521 247.780829 \n",
       "L 361.385346 249.719052 \n",
       "L 362.444037 250.725659 \n",
       "L 363.522121 251.114749 \n",
       "L 364.61286 251.154957 \n",
       "L 365.711262 251.005768 \n",
       "L 367.924204 250.383627 \n",
       "L 371.266786 249.107424 \n",
       "L 373.519987 248.004832 \n",
       "L 374.659467 247.342975 \n",
       "L 376.989941 245.602023 \n",
       "L 378.193256 244.44484 \n",
       "L 379.430647 243.050076 \n",
       "L 382.005868 239.661794 \n",
       "L 383.3152 237.930135 \n",
       "L 384.599111 236.476734 \n",
       "L 385.815212 235.584017 \n",
       "L 386.931877 235.424019 \n",
       "L 387.95627 235.88255 \n",
       "L 389.830703 237.785399 \n",
       "L 390.757202 238.694883 \n",
       "L 391.701899 239.42409 \n",
       "L 392.681555 239.886671 \n",
       "L 393.687234 240.154793 \n",
       "L 394.719415 240.23991 \n",
       "L 396.840466 240.035844 \n",
       "L 400.103083 239.22344 \n",
       "L 407.827518 236.672837 \n",
       "L 471.90612 214.401936 \n",
       "L 492.842176 207.115891 \n",
       "L 508.241227 201.770718 \n",
       "L 535.716498 192.208963 \n",
       "L 551.070241 186.878549 \n",
       "L 558.74845 184.204241 \n",
       "L 558.74845 184.204241 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 268.715724 276.67643 \n",
       "L 342.234069 251.021895 \n",
       "L 344.450017 249.981141 \n",
       "L 345.555589 249.278909 \n",
       "L 346.659275 248.315866 \n",
       "L 347.762117 246.901149 \n",
       "L 348.867621 244.725514 \n",
       "L 349.984554 241.381338 \n",
       "L 351.129103 236.564198 \n",
       "L 352.319776 230.61477 \n",
       "L 353.554757 225.342116 \n",
       "L 354.786043 223.481952 \n",
       "L 355.947355 226.079082 \n",
       "L 358.057176 235.985555 \n",
       "L 359.086091 239.485162 \n",
       "L 360.131831 241.56805 \n",
       "L 361.196715 242.644534 \n",
       "L 362.277708 243.078963 \n",
       "L 363.371089 243.118083 \n",
       "L 364.472098 242.946378 \n",
       "L 366.694039 242.194779 \n",
       "L 368.941903 241.074194 \n",
       "L 370.080069 240.356046 \n",
       "L 371.233677 239.481581 \n",
       "L 372.408005 238.41255 \n",
       "L 373.607473 237.12756 \n",
       "L 376.099579 233.841346 \n",
       "L 378.655506 230.355681 \n",
       "L 379.89047 229.154704 \n",
       "L 381.055152 228.622473 \n",
       "L 382.133319 228.823593 \n",
       "L 383.138123 229.576584 \n",
       "L 386.009013 232.460647 \n",
       "L 386.990626 233.084027 \n",
       "L 387.998734 233.465234 \n",
       "L 389.031314 233.638383 \n",
       "L 390.082114 233.661446 \n",
       "L 392.22848 233.365293 \n",
       "L 394.408797 232.820607 \n",
       "L 397.707609 231.806193 \n",
       "L 406.550941 228.790172 \n",
       "L 550.790118 178.682679 \n",
       "L 550.790118 178.682679 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 260.844475 271.051077 \n",
       "L 287.616321 261.764656 \n",
       "L 303.20792 256.343408 \n",
       "L 331.014064 246.696428 \n",
       "L 341.006234 243.103919 \n",
       "L 343.219534 242.064921 \n",
       "L 344.322773 241.357111 \n",
       "L 345.42318 240.398946 \n",
       "L 346.520166 238.971191 \n",
       "L 347.616212 236.792083 \n",
       "L 348.717485 233.434107 \n",
       "L 349.838169 228.563433 \n",
       "L 350.998426 222.512339 \n",
       "L 352.206466 217.004158 \n",
       "L 353.430242 214.872062 \n",
       "L 354.609067 217.319793 \n",
       "L 356.775872 227.466658 \n",
       "L 357.824232 231.126052 \n",
       "L 358.881839 233.324328 \n",
       "L 359.954473 234.445839 \n",
       "L 361.041816 234.85744 \n",
       "L 362.141167 234.833414 \n",
       "L 363.251323 234.516935 \n",
       "L 364.372269 233.980703 \n",
       "L 365.505281 233.252128 \n",
       "L 366.654092 232.310183 \n",
       "L 367.822704 231.134966 \n",
       "L 369.015762 229.701924 \n",
       "L 371.478682 226.184246 \n",
       "L 372.73125 224.376674 \n",
       "L 373.967859 222.891039 \n",
       "L 375.159217 222.006658 \n",
       "L 376.283848 221.871066 \n",
       "L 377.339446 222.412928 \n",
       "L 379.324596 224.456217 \n",
       "L 380.303735 225.440138 \n",
       "L 381.297171 226.200934 \n",
       "L 382.309935 226.719777 \n",
       "L 383.34403 227.004743 \n",
       "L 384.395814 227.107128 \n",
       "L 386.540081 226.921741 \n",
       "L 388.717675 226.435935 \n",
       "L 393.118085 225.089292 \n",
       "L 403.065866 221.68275 \n",
       "L 486.872906 192.586822 \n",
       "L 496.75943 189.166245 \n",
       "L 507.743894 185.352698 \n",
       "L 524.209024 179.631196 \n",
       "L 542.843929 173.161464 \n",
       "L 542.843929 173.161464 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 252.981766 265.431976 \n",
       "L 339.779884 235.189928 \n",
       "L 341.990229 234.139939 \n",
       "L 343.090678 233.407742 \n",
       "L 344.187721 232.43724 \n",
       "L 345.278984 230.982805 \n",
       "L 346.365846 228.770154 \n",
       "L 347.452631 225.393762 \n",
       "L 348.551075 220.495062 \n",
       "L 349.682552 214.42578 \n",
       "L 350.863118 208.842195 \n",
       "L 352.076159 206.530351 \n",
       "L 353.269138 208.789054 \n",
       "L 355.495102 218.848508 \n",
       "L 356.566318 222.537991 \n",
       "L 357.640539 224.714501 \n",
       "L 358.726526 225.73989 \n",
       "L 359.828321 225.931493 \n",
       "L 360.945859 225.571915 \n",
       "L 362.081517 224.770422 \n",
       "L 363.238247 223.583154 \n",
       "L 364.417031 222.086279 \n",
       "L 368.059455 216.768171 \n",
       "L 369.25983 215.54085 \n",
       "L 370.414352 215.057613 \n",
       "L 371.5145 215.310164 \n",
       "L 372.563866 216.138857 \n",
       "L 374.586468 218.326399 \n",
       "L 375.595079 219.250182 \n",
       "L 376.618477 219.905981 \n",
       "L 377.657806 220.315294 \n",
       "L 378.711004 220.527229 \n",
       "L 380.85664 220.466305 \n",
       "L 383.03492 220.033371 \n",
       "L 385.22939 219.42839 \n",
       "L 389.63889 218.01364 \n",
       "L 441.520048 200.048996 \n",
       "L 451.441038 196.592305 \n",
       "L 461.343477 193.172732 \n",
       "L 471.25199 189.72136 \n",
       "L 482.244504 185.913029 \n",
       "L 534.898027 167.660247 \n",
       "L 534.898027 167.660247 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 245.130201 259.824755 \n",
       "L 337.447702 227.708994 \n",
       "L 339.656677 226.717404 \n",
       "L 340.758294 226.114816 \n",
       "L 341.855042 225.321443 \n",
       "L 342.94686 224.255811 \n",
       "L 344.031753 222.710911 \n",
       "L 345.108861 220.373369 \n",
       "L 346.182654 216.933885 \n",
       "L 347.263127 212.110401 \n",
       "L 348.371703 206.336762 \n",
       "L 349.528692 201.231627 \n",
       "L 350.728113 199.136522 \n",
       "L 351.929029 201.094752 \n",
       "L 354.215928 209.912081 \n",
       "L 355.318201 213.115208 \n",
       "L 356.420437 214.801739 \n",
       "L 357.53436 215.269789 \n",
       "L 358.667901 214.774748 \n",
       "L 359.823485 213.611221 \n",
       "L 363.377134 208.987216 \n",
       "L 364.544393 208.187716 \n",
       "L 365.675036 208.15216 \n",
       "L 366.762516 208.833289 \n",
       "L 367.818112 209.900505 \n",
       "L 369.884511 212.176476 \n",
       "L 370.922462 212.9938 \n",
       "L 371.968889 213.57281 \n",
       "L 373.029404 213.878995 \n",
       "L 374.099169 214.008855 \n",
       "L 375.17777 213.990996 \n",
       "L 377.355462 213.648247 \n",
       "L 380.649746 212.754003 \n",
       "L 395.004247 207.888499 \n",
       "L 402.735352 205.204205 \n",
       "L 526.976318 162.137231 \n",
       "L 526.976318 162.137231 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 237.283696 254.216055 \n",
       "L 310.724399 228.782073 \n",
       "L 326.256689 223.404965 \n",
       "L 335.113467 220.141152 \n",
       "L 337.317801 219.108008 \n",
       "L 338.414972 218.451472 \n",
       "L 339.508889 217.667011 \n",
       "L 340.59713 216.650161 \n",
       "L 341.67771 215.258973 \n",
       "L 342.75131 213.362748 \n",
       "L 343.817751 210.715801 \n",
       "L 344.88371 207.203217 \n",
       "L 347.070504 198.618252 \n",
       "L 348.222187 195.45542 \n",
       "L 349.408381 194.493124 \n",
       "L 350.602334 195.761761 \n",
       "L 352.936621 200.246582 \n",
       "L 354.082444 201.471873 \n",
       "L 355.230758 201.624685 \n",
       "L 356.388656 201.017618 \n",
       "L 357.554252 200.138042 \n",
       "L 358.717476 199.488519 \n",
       "L 359.861548 199.530369 \n",
       "L 360.976939 200.302734 \n",
       "L 362.061468 201.637922 \n",
       "L 364.175888 204.604535 \n",
       "L 365.226709 205.773702 \n",
       "L 366.284147 206.59365 \n",
       "L 367.349863 207.106802 \n",
       "L 368.422896 207.387091 \n",
       "L 369.50509 207.44793 \n",
       "L 371.683929 207.228583 \n",
       "L 374.978904 206.391271 \n",
       "L 380.493562 204.605842 \n",
       "L 519.047227 156.649416 \n",
       "L 519.047227 156.649416 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 229.447436 248.617135 \n",
       "L 259.528312 238.213811 \n",
       "L 271.762568 233.97102 \n",
       "L 305.088126 222.451344 \n",
       "L 315.068308 218.98215 \n",
       "L 325.04284 215.51389 \n",
       "L 330.568867 213.4127 \n",
       "L 332.768186 212.381056 \n",
       "L 333.86435 211.784619 \n",
       "L 334.954766 211.060316 \n",
       "L 336.040488 210.198839 \n",
       "L 337.119827 209.132877 \n",
       "L 338.190902 207.767258 \n",
       "L 339.256625 206.085368 \n",
       "L 340.317227 203.99273 \n",
       "L 343.552357 196.622873 \n",
       "L 344.678228 194.849612 \n",
       "L 345.829167 193.88273 \n",
       "L 346.993995 193.608506 \n",
       "L 348.160197 193.587414 \n",
       "L 349.321144 193.194644 \n",
       "L 350.480528 192.066405 \n",
       "L 351.647266 190.062178 \n",
       "L 352.828947 187.629547 \n",
       "L 354.021309 185.709178 \n",
       "L 355.204888 185.283234 \n",
       "L 356.354367 186.802959 \n",
       "L 357.460126 189.665337 \n",
       "L 359.585607 195.747078 \n",
       "L 360.639543 197.85046 \n",
       "L 361.69892 199.284686 \n",
       "L 362.768296 200.12776 \n",
       "L 363.846329 200.563062 \n",
       "L 364.930832 200.733109 \n",
       "L 367.116355 200.546763 \n",
       "L 369.314145 200.023156 \n",
       "L 372.619522 199.021747 \n",
       "L 379.240589 196.789446 \n",
       "L 389.174627 193.342848 \n",
       "L 399.098897 189.921928 \n",
       "L 410.120337 186.106969 \n",
       "L 432.14109 178.477581 \n",
       "L 442.035287 175.065143 \n",
       "L 453.024649 171.262647 \n",
       "L 468.399757 165.933804 \n",
       "L 511.138097 151.145011 \n",
       "L 511.138097 151.145011 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 221.620711 243.026218 \n",
       "L 234.992254 238.401697 \n",
       "L 321.615502 208.398621 \n",
       "L 326.025981 206.648986 \n",
       "L 328.218463 205.58923 \n",
       "L 330.39026 204.183837 \n",
       "L 331.463829 203.247157 \n",
       "L 332.529224 202.116908 \n",
       "L 333.584907 200.729899 \n",
       "L 334.630719 199.02984 \n",
       "L 337.789106 193.023526 \n",
       "L 338.886744 191.504695 \n",
       "L 340.018777 190.681544 \n",
       "L 341.181777 190.697338 \n",
       "L 342.358352 191.291395 \n",
       "L 343.532977 192.11519 \n",
       "L 344.693928 192.744524 \n",
       "L 345.836032 192.755104 \n",
       "L 346.962316 191.846128 \n",
       "L 348.079474 189.625097 \n",
       "L 349.201027 185.845394 \n",
       "L 351.525618 175.353939 \n",
       "L 352.731901 172.177854 \n",
       "L 353.923152 172.944546 \n",
       "L 355.06153 177.122828 \n",
       "L 356.147214 182.38875 \n",
       "L 357.20597 186.88352 \n",
       "L 358.260469 190.032678 \n",
       "L 359.323209 191.93339 \n",
       "L 360.396322 192.96292 \n",
       "L 361.478411 193.437713 \n",
       "L 362.568374 193.54516 \n",
       "L 363.662229 193.474511 \n",
       "L 365.860202 192.985527 \n",
       "L 369.166115 191.971444 \n",
       "L 379.095384 188.577179 \n",
       "L 414.355192 176.382895 \n",
       "L 424.25401 172.971844 \n",
       "L 434.152074 169.537903 \n",
       "L 447.331979 164.991094 \n",
       "L 460.504255 160.433707 \n",
       "L 503.229489 145.663109 \n",
       "L 503.229489 145.663109 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 213.796008 237.429333 \n",
       "L 220.489734 235.134755 \n",
       "L 228.278391 232.426134 \n",
       "L 234.964834 230.134193 \n",
       "L 242.747724 227.43007 \n",
       "L 250.531195 224.737156 \n",
       "L 263.873374 220.148594 \n",
       "L 270.532157 217.831979 \n",
       "L 319.291092 200.834864 \n",
       "L 322.579351 199.349289 \n",
       "L 323.66824 198.760679 \n",
       "L 324.74758 198.049369 \n",
       "L 325.817064 197.195148 \n",
       "L 326.877093 196.182666 \n",
       "L 327.920274 194.894411 \n",
       "L 328.952031 193.360149 \n",
       "L 330.988685 189.494843 \n",
       "L 332.019915 187.484795 \n",
       "L 333.078766 185.731963 \n",
       "L 334.184954 184.646846 \n",
       "L 335.339159 184.39241 \n",
       "L 336.529009 184.909493 \n",
       "L 338.937497 187.188036 \n",
       "L 340.122834 188.290817 \n",
       "L 341.286479 189.076695 \n",
       "L 342.428073 189.406875 \n",
       "L 343.552175 189.256581 \n",
       "L 344.659345 188.396492 \n",
       "L 345.754034 186.625193 \n",
       "L 346.840976 183.546304 \n",
       "L 347.930342 178.722646 \n",
       "L 350.19848 165.670069 \n",
       "L 351.397479 161.822902 \n",
       "L 352.597476 163.074224 \n",
       "L 353.751511 168.31571 \n",
       "L 354.852564 174.503057 \n",
       "L 355.924769 179.454902 \n",
       "L 356.990752 182.710878 \n",
       "L 358.061954 184.584244 \n",
       "L 359.1419 185.510945 \n",
       "L 360.228804 185.89154 \n",
       "L 361.321561 185.934212 \n",
       "L 363.516962 185.549653 \n",
       "L 367.924475 184.16773 \n",
       "L 375.642781 181.525077 \n",
       "L 408.680936 170.115518 \n",
       "L 495.333129 140.17894 \n",
       "L 495.333129 140.17894 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 205.992664 231.861997 \n",
       "L 210.44211 230.311506 \n",
       "L 217.133101 228.022464 \n",
       "L 221.572778 226.458972 \n",
       "L 226.0363 224.946932 \n",
       "L 230.476975 223.390025 \n",
       "L 234.932725 221.869482 \n",
       "L 241.595109 219.545983 \n",
       "L 246.051451 218.038243 \n",
       "L 250.486381 216.479178 \n",
       "L 254.936916 214.966256 \n",
       "L 259.373364 213.417648 \n",
       "L 266.03976 211.138664 \n",
       "L 270.470321 209.580117 \n",
       "L 274.914187 208.072621 \n",
       "L 279.344859 206.522038 \n",
       "L 283.782787 205.005336 \n",
       "L 290.428199 202.692721 \n",
       "L 297.078152 200.416659 \n",
       "L 314.764674 194.124847 \n",
       "L 318.038268 192.561588 \n",
       "L 320.185125 191.128355 \n",
       "L 321.237676 190.164858 \n",
       "L 322.275897 189.006218 \n",
       "L 323.296412 187.591547 \n",
       "L 325.287962 183.899323 \n",
       "L 327.295984 179.848333 \n",
       "L 328.361784 178.416709 \n",
       "L 329.486473 177.72383 \n",
       "L 330.673509 177.950442 \n",
       "L 331.896643 178.839078 \n",
       "L 335.557758 182.489907 \n",
       "L 336.736262 183.327101 \n",
       "L 337.891494 183.836719 \n",
       "L 339.027983 184.040619 \n",
       "L 340.151708 184.010831 \n",
       "L 341.261808 183.639486 \n",
       "L 342.362328 182.936378 \n",
       "L 343.451944 181.700488 \n",
       "L 344.529507 179.600761 \n",
       "L 345.598078 176.272918 \n",
       "L 346.663156 171.149446 \n",
       "L 348.871054 157.137675 \n",
       "L 350.054238 152.817387 \n",
       "L 351.261896 154.022871 \n",
       "L 352.43948 159.578332 \n",
       "L 353.566401 166.146264 \n",
       "L 354.658232 171.416736 \n",
       "L 355.737593 174.814096 \n",
       "L 356.816499 176.775441 \n",
       "L 357.900778 177.733356 \n",
       "L 358.99049 178.104108 \n",
       "L 360.083846 178.163993 \n",
       "L 362.280151 177.753357 \n",
       "L 365.582759 176.75148 \n",
       "L 382.108879 171.058098 \n",
       "L 392.01511 167.645773 \n",
       "L 400.816916 164.600351 \n",
       "L 441.455879 150.563292 \n",
       "L 445.834988 149.075953 \n",
       "L 450.226557 147.539531 \n",
       "L 456.794611 145.293789 \n",
       "L 463.372241 143.008059 \n",
       "L 487.444732 134.701726 \n",
       "L 487.444732 134.701726 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 198.181568 226.271965 \n",
       "L 202.640697 224.742478 \n",
       "L 212.65758 221.282904 \n",
       "L 221.56129 218.220452 \n",
       "L 305.840933 189.103857 \n",
       "L 311.331604 186.902335 \n",
       "L 313.495278 185.737758 \n",
       "L 314.561863 185.010315 \n",
       "L 315.613239 184.130132 \n",
       "L 316.650028 183.087407 \n",
       "L 317.659782 181.750718 \n",
       "L 318.650402 180.163938 \n",
       "L 320.578176 176.166124 \n",
       "L 321.553657 174.099483 \n",
       "L 322.567616 172.322958 \n",
       "L 323.652049 171.245951 \n",
       "L 324.81412 171.050107 \n",
       "L 326.037422 171.651432 \n",
       "L 328.556418 174.179803 \n",
       "L 329.797777 175.447439 \n",
       "L 331.01128 176.47094 \n",
       "L 332.194881 177.171412 \n",
       "L 333.356648 177.628138 \n",
       "L 334.498466 177.826046 \n",
       "L 335.626891 177.840221 \n",
       "L 336.746267 177.721593 \n",
       "L 338.961851 177.039245 \n",
       "L 340.060104 176.442518 \n",
       "L 341.150418 175.54802 \n",
       "L 342.231847 174.208373 \n",
       "L 343.299734 172.038946 \n",
       "L 344.354278 168.62742 \n",
       "L 345.399466 163.445564 \n",
       "L 347.54847 148.944548 \n",
       "L 348.712442 144.250532 \n",
       "L 349.92479 145.199375 \n",
       "L 351.126416 150.853277 \n",
       "L 352.281376 157.646737 \n",
       "L 353.395156 163.186574 \n",
       "L 354.488426 166.776442 \n",
       "L 355.575737 168.830952 \n",
       "L 356.663763 169.881101 \n",
       "L 357.755805 170.266452 \n",
       "L 358.849953 170.346813 \n",
       "L 359.947131 170.202162 \n",
       "L 362.145236 169.658643 \n",
       "L 379.759986 163.657008 \n",
       "L 388.561686 160.614134 \n",
       "L 414.935792 151.506866 \n",
       "L 419.320723 150.022191 \n",
       "L 423.715319 148.483221 \n",
       "L 430.291487 146.237965 \n",
       "L 437.97293 143.563972 \n",
       "L 442.351894 142.071337 \n",
       "L 448.93157 139.781504 \n",
       "L 454.406791 137.895158 \n",
       "L 462.061491 135.278108 \n",
       "L 466.44496 133.7476 \n",
       "L 472.99864 131.511205 \n",
       "L 477.383801 129.969251 \n",
       "L 479.566132 129.226544 \n",
       "L 479.566132 129.226544 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 190.392109 220.710863 \n",
       "L 202.631559 216.481497 \n",
       "L 212.65003 213.045779 \n",
       "L 222.648157 209.585287 \n",
       "L 235.981636 205.003093 \n",
       "L 245.965172 201.547255 \n",
       "L 254.839191 198.490891 \n",
       "L 264.818477 195.063001 \n",
       "L 274.784899 191.616149 \n",
       "L 284.750287 188.193456 \n",
       "L 304.625857 181.13593 \n",
       "L 306.803122 180.134774 \n",
       "L 308.948738 178.876463 \n",
       "L 310.00274 178.08471 \n",
       "L 311.035859 177.106947 \n",
       "L 312.044724 175.899531 \n",
       "L 313.027809 174.428845 \n",
       "L 313.981076 172.635661 \n",
       "L 316.81524 166.433573 \n",
       "L 317.846974 164.977402 \n",
       "L 318.965664 164.324023 \n",
       "L 320.173492 164.588838 \n",
       "L 321.441693 165.566859 \n",
       "L 325.263166 169.454132 \n",
       "L 326.48142 170.345956 \n",
       "L 327.668996 170.946121 \n",
       "L 328.830477 171.27492 \n",
       "L 329.97439 171.416701 \n",
       "L 331.103967 171.390852 \n",
       "L 333.338568 171.038577 \n",
       "L 335.553047 170.398241 \n",
       "L 337.75328 169.46414 \n",
       "L 338.847119 168.81162 \n",
       "L 339.933031 167.87748 \n",
       "L 341.008847 166.497671 \n",
       "L 342.069959 164.337997 \n",
       "L 343.112502 160.904186 \n",
       "L 344.13936 155.700046 \n",
       "L 346.231001 140.857054 \n",
       "L 347.37385 135.821583 \n",
       "L 348.587909 136.467271 \n",
       "L 349.81245 142.116127 \n",
       "L 350.995803 149.124295 \n",
       "L 352.132691 154.892332 \n",
       "L 353.240674 158.701616 \n",
       "L 354.336615 160.867887 \n",
       "L 355.429131 161.997283 \n",
       "L 356.522932 162.437279 \n",
       "L 357.618238 162.524042 \n",
       "L 358.715135 162.406033 \n",
       "L 360.912786 161.859832 \n",
       "L 366.414062 160.046384 \n",
       "L 471.691848 123.769014 \n",
       "L 471.691848 123.769014 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 182.60346 215.14295 \n",
       "L 203.744835 207.870546 \n",
       "L 219.298327 202.508603 \n",
       "L 296.824564 175.734705 \n",
       "L 301.198883 173.92589 \n",
       "L 303.345937 172.737796 \n",
       "L 304.400997 172.003075 \n",
       "L 305.435235 171.105672 \n",
       "L 306.444388 170.002389 \n",
       "L 307.425336 168.653842 \n",
       "L 308.36976 166.97715 \n",
       "L 309.288737 165.030757 \n",
       "L 311.10772 160.705534 \n",
       "L 312.080626 158.918831 \n",
       "L 313.143094 157.828119 \n",
       "L 314.313863 157.670416 \n",
       "L 315.573797 158.360812 \n",
       "L 316.882157 159.593845 \n",
       "L 319.491921 162.349716 \n",
       "L 320.74967 163.424378 \n",
       "L 321.968587 164.181398 \n",
       "L 323.154422 164.648174 \n",
       "L 324.316275 164.897573 \n",
       "L 325.457993 164.949081 \n",
       "L 327.70793 164.716139 \n",
       "L 331.037584 163.86524 \n",
       "L 334.34601 162.731187 \n",
       "L 336.541129 161.748901 \n",
       "L 337.631633 161.060278 \n",
       "L 338.715639 160.144369 \n",
       "L 339.786808 158.747315 \n",
       "L 340.841669 156.596968 \n",
       "L 341.873994 153.18598 \n",
       "L 342.882943 147.947555 \n",
       "L 344.918684 132.830722 \n",
       "L 346.038487 127.438616 \n",
       "L 347.251889 127.791939 \n",
       "L 348.497536 133.370168 \n",
       "L 349.709528 140.593938 \n",
       "L 350.870114 146.589091 \n",
       "L 351.993787 150.602797 \n",
       "L 353.098617 152.917788 \n",
       "L 354.196172 154.102226 \n",
       "L 355.291832 154.616596 \n",
       "L 356.388211 154.712672 \n",
       "L 357.485205 154.610438 \n",
       "L 359.682187 154.069662 \n",
       "L 362.980386 153.015473 \n",
       "L 371.776627 150.003036 \n",
       "L 392.647422 142.823028 \n",
       "L 402.52557 139.408589 \n",
       "L 411.298373 136.392245 \n",
       "L 426.638123 131.116305 \n",
       "L 439.776122 126.587494 \n",
       "L 463.829655 118.306575 \n",
       "L 463.829655 118.306575 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 174.830533 209.592878 \n",
       "L 185.954693 205.762424 \n",
       "L 237.021103 188.175782 \n",
       "L 243.676789 185.906915 \n",
       "L 253.638634 182.460124 \n",
       "L 270.240948 176.764227 \n",
       "L 276.870329 174.462649 \n",
       "L 282.396909 172.562462 \n",
       "L 290.126155 169.880061 \n",
       "L 295.602854 167.719006 \n",
       "L 297.752027 166.590938 \n",
       "L 298.808984 165.906253 \n",
       "L 299.843233 165.064479 \n",
       "L 300.859557 164.08566 \n",
       "L 301.837034 162.821316 \n",
       "L 302.781914 161.293148 \n",
       "L 303.690483 159.455761 \n",
       "L 306.358828 153.049384 \n",
       "L 307.356719 151.561727 \n",
       "L 308.473897 150.96091 \n",
       "L 309.701174 151.248689 \n",
       "L 311.017438 152.326344 \n",
       "L 314.996997 156.41805 \n",
       "L 316.250809 157.32336 \n",
       "L 317.469275 157.969313 \n",
       "L 318.652194 158.325968 \n",
       "L 320.951868 158.483294 \n",
       "L 323.197074 158.147564 \n",
       "L 326.52647 157.279796 \n",
       "L 333.137066 154.994536 \n",
       "L 335.329377 154.005621 \n",
       "L 336.417236 153.297047 \n",
       "L 337.499015 152.384852 \n",
       "L 338.5673 151.013424 \n",
       "L 339.615265 148.84399 \n",
       "L 340.638664 145.475717 \n",
       "L 341.630794 140.222157 \n",
       "L 343.61179 124.864277 \n",
       "L 344.706327 119.084444 \n",
       "L 345.916892 119.145261 \n",
       "L 347.181883 124.637634 \n",
       "L 348.4224 132.040318 \n",
       "L 349.607426 138.298673 \n",
       "L 350.747456 142.493168 \n",
       "L 351.861828 144.969516 \n",
       "L 352.964463 146.228715 \n",
       "L 354.062674 146.781481 \n",
       "L 355.159651 146.931449 \n",
       "L 357.35463 146.610372 \n",
       "L 360.649993 145.630744 \n",
       "L 384.813623 137.326672 \n",
       "L 391.395248 135.07664 \n",
       "L 425.362298 123.385215 \n",
       "L 451.603771 114.36277 \n",
       "L 455.976804 112.846046 \n",
       "L 455.976804 112.846046 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 167.04846 204.02128 \n",
       "L 171.515844 202.51856 \n",
       "L 175.94853 200.963044 \n",
       "L 180.409745 199.457019 \n",
       "L 184.847004 197.913266 \n",
       "L 191.524658 195.64119 \n",
       "L 195.953123 194.086968 \n",
       "L 200.407604 192.585151 \n",
       "L 204.837061 191.037247 \n",
       "L 209.284123 189.528225 \n",
       "L 213.715167 187.9887 \n",
       "L 220.377194 185.719846 \n",
       "L 224.799265 184.166495 \n",
       "L 229.243068 182.669189 \n",
       "L 233.666051 181.123447 \n",
       "L 238.102257 179.616028 \n",
       "L 242.528187 178.084815 \n",
       "L 249.174531 175.817767 \n",
       "L 253.592001 174.270244 \n",
       "L 258.023214 172.771635 \n",
       "L 262.441528 171.234373 \n",
       "L 266.866185 169.72469 \n",
       "L 273.491461 167.425077 \n",
       "L 277.916844 165.935311 \n",
       "L 282.324813 164.373946 \n",
       "L 285.632898 163.213932 \n",
       "L 288.924315 161.967829 \n",
       "L 292.166901 160.440565 \n",
       "L 294.266532 159.037068 \n",
       "L 295.280404 158.102372 \n",
       "L 296.263416 156.957278 \n",
       "L 297.214222 155.579724 \n",
       "L 298.112073 153.822811 \n",
       "L 298.981296 151.818414 \n",
       "L 300.684763 147.342994 \n",
       "L 301.619813 145.550683 \n",
       "L 302.659548 144.435309 \n",
       "L 303.843844 144.339929 \n",
       "L 305.1396 145.096063 \n",
       "L 306.497763 146.394163 \n",
       "L 309.215173 149.285038 \n",
       "L 310.516447 150.409499 \n",
       "L 311.770111 151.214012 \n",
       "L 312.977439 151.680488 \n",
       "L 314.160233 151.976366 \n",
       "L 316.449458 151.960815 \n",
       "L 317.576923 151.819626 \n",
       "L 320.913034 150.990444 \n",
       "L 328.632091 148.473205 \n",
       "L 333.023832 146.758092 \n",
       "L 334.118383 146.248256 \n",
       "L 335.205245 145.554864 \n",
       "L 336.283893 144.623924 \n",
       "L 337.349804 143.282909 \n",
       "L 338.391924 141.118541 \n",
       "L 339.405701 137.753075 \n",
       "L 340.383131 132.531414 \n",
       "L 342.309786 116.924015 \n",
       "L 343.378063 110.789102 \n",
       "L 344.582704 110.50101 \n",
       "L 345.865901 115.93424 \n",
       "L 348.344616 130.000815 \n",
       "L 349.501629 134.395559 \n",
       "L 350.626167 137.007398 \n",
       "L 351.734146 138.376523 \n",
       "L 352.835052 138.955504 \n",
       "L 353.933025 139.1497 \n",
       "L 355.030666 139.050141 \n",
       "L 357.225516 138.565258 \n",
       "L 368.205042 134.89712 \n",
       "L 372.597148 133.357885 \n",
       "L 376.984349 131.878694 \n",
       "L 381.37438 130.342711 \n",
       "L 387.951625 128.104059 \n",
       "L 395.626692 125.441931 \n",
       "L 400.006581 123.955109 \n",
       "L 406.580797 121.675053 \n",
       "L 413.146808 119.431748 \n",
       "L 433.925106 112.296724 \n",
       "L 438.300943 110.772979 \n",
       "L 442.667139 109.286155 \n",
       "L 448.125346 107.417514 \n",
       "L 448.125346 107.417514 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 159.306551 198.50474 \n",
       "L 163.737824 196.9495 \n",
       "L 168.200191 195.445183 \n",
       "L 172.63503 193.899371 \n",
       "L 177.086374 192.383284 \n",
       "L 183.741738 190.077018 \n",
       "L 188.197128 188.57632 \n",
       "L 192.624358 187.026556 \n",
       "L 197.073745 185.52187 \n",
       "L 201.503154 183.980964 \n",
       "L 205.944332 182.466989 \n",
       "L 215.918174 179.032874 \n",
       "L 222.573481 176.767829 \n",
       "L 226.992402 175.218351 \n",
       "L 231.430976 173.722287 \n",
       "L 235.850471 172.179988 \n",
       "L 240.281416 170.672911 \n",
       "L 246.911294 168.372165 \n",
       "L 251.344011 166.883842 \n",
       "L 255.756717 165.33561 \n",
       "L 259.077149 164.214249 \n",
       "L 271.229725 160.052292 \n",
       "L 275.636412 158.49695 \n",
       "L 280.049467 156.972803 \n",
       "L 285.51197 154.801495 \n",
       "L 287.646648 153.65549 \n",
       "L 288.692314 152.955348 \n",
       "L 289.71009 152.0901 \n",
       "L 290.706137 151.084628 \n",
       "L 291.653164 149.774681 \n",
       "L 292.560321 148.195215 \n",
       "L 293.421335 146.292374 \n",
       "L 295.926156 139.692766 \n",
       "L 296.893581 138.191269 \n",
       "L 298.008628 137.618449 \n",
       "L 299.260006 137.96243 \n",
       "L 300.623076 139.117266 \n",
       "L 304.758921 143.404251 \n",
       "L 306.049284 144.346543 \n",
       "L 307.297036 145.024605 \n",
       "L 308.499909 145.39834 \n",
       "L 309.671126 145.551932 \n",
       "L 310.823893 145.579156 \n",
       "L 313.079639 145.258659 \n",
       "L 316.41406 144.402821 \n",
       "L 327.427319 140.717523 \n",
       "L 330.72241 139.507553 \n",
       "L 332.908505 138.487572 \n",
       "L 333.995351 137.826289 \n",
       "L 335.070974 136.879604 \n",
       "L 336.134055 135.551908 \n",
       "L 337.171443 133.413463 \n",
       "L 338.175766 130.044554 \n",
       "L 339.139062 124.84643 \n",
       "L 340.068109 117.465732 \n",
       "L 341.012276 108.997221 \n",
       "L 342.053779 102.541477 \n",
       "L 343.249638 101.876266 \n",
       "L 344.549446 107.232487 \n",
       "L 347.081524 121.686903 \n",
       "L 348.256412 126.312183 \n",
       "L 349.391483 129.046408 \n",
       "L 350.505315 130.525578 \n",
       "L 351.608939 131.15029 \n",
       "L 352.708217 131.363616 \n",
       "L 353.805896 131.302808 \n",
       "L 356.000149 130.822767 \n",
       "L 364.778238 127.914824 \n",
       "L 369.167416 126.380852 \n",
       "L 373.552953 124.895539 \n",
       "L 380.132152 122.614558 \n",
       "L 387.80007 120.005671 \n",
       "L 392.183253 118.481602 \n",
       "L 398.749696 116.248846 \n",
       "L 403.132308 114.710253 \n",
       "L 407.504359 113.238968 \n",
       "L 411.884053 111.708295 \n",
       "L 416.255106 110.224299 \n",
       "L 422.818717 107.944941 \n",
       "L 427.182712 106.479231 \n",
       "L 431.559338 104.943912 \n",
       "L 435.921017 103.474392 \n",
       "L 440.293924 101.948198 \n",
       "L 440.293924 101.948198 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 151.537616 192.940429 \n",
       "L 157.098332 191.034071 \n",
       "L 164.885242 188.375277 \n",
       "L 174.876658 184.93584 \n",
       "L 184.873516 181.519931 \n",
       "L 191.520738 179.218492 \n",
       "L 202.610974 175.419355 \n",
       "L 220.329755 169.328439 \n",
       "L 230.295528 165.926223 \n",
       "L 240.242384 162.493403 \n",
       "L 251.29621 158.709077 \n",
       "L 278.850276 149.062923 \n",
       "L 282.076742 147.517006 \n",
       "L 284.153486 146.083723 \n",
       "L 285.150606 145.134219 \n",
       "L 286.108332 143.954983 \n",
       "L 287.023512 142.516265 \n",
       "L 287.880155 140.715292 \n",
       "L 289.482617 136.313728 \n",
       "L 290.288549 134.033187 \n",
       "L 291.179929 132.186364 \n",
       "L 292.206455 131.102649 \n",
       "L 293.401269 131.03605 \n",
       "L 294.73493 131.864885 \n",
       "L 296.146244 133.250755 \n",
       "L 298.969395 136.269063 \n",
       "L 300.312406 137.435163 \n",
       "L 301.59727 138.262338 \n",
       "L 302.830803 138.775377 \n",
       "L 304.028057 139.062952 \n",
       "L 305.192929 139.137623 \n",
       "L 307.469918 138.937028 \n",
       "L 310.812382 138.123815 \n",
       "L 316.328064 136.366872 \n",
       "L 329.516336 131.755174 \n",
       "L 331.700935 130.747639 \n",
       "L 332.786368 130.087982 \n",
       "L 333.860999 129.166988 \n",
       "L 334.919884 127.819387 \n",
       "L 335.953506 125.720733 \n",
       "L 336.949129 122.358736 \n",
       "L 337.898193 117.161139 \n",
       "L 338.805819 109.76255 \n",
       "L 339.719424 101.09514 \n",
       "L 340.733299 94.326397 \n",
       "L 341.918322 93.292899 \n",
       "L 343.232331 98.520656 \n",
       "L 345.818028 113.358466 \n",
       "L 347.01174 118.227569 \n",
       "L 348.157834 121.100638 \n",
       "L 349.277839 122.666592 \n",
       "L 350.384447 123.367792 \n",
       "L 351.485047 123.586341 \n",
       "L 352.583009 123.560466 \n",
       "L 354.776772 123.077188 \n",
       "L 358.06644 122.055624 \n",
       "L 367.934602 118.663497 \n",
       "L 375.603977 116.049774 \n",
       "L 382.177311 113.770206 \n",
       "L 388.744885 111.528361 \n",
       "L 397.497574 108.531701 \n",
       "L 432.459445 96.533611 \n",
       "L 432.459445 96.533611 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 143.804112 187.422489 \n",
       "L 157.131396 182.837066 \n",
       "L 167.132677 179.422204 \n",
       "L 177.112122 175.984075 \n",
       "L 187.099329 172.574203 \n",
       "L 194.851509 169.902444 \n",
       "L 207.038817 165.734107 \n",
       "L 218.102514 161.932115 \n",
       "L 265.592807 145.632207 \n",
       "L 273.278505 142.844206 \n",
       "L 276.511814 141.363222 \n",
       "L 278.599862 140.033825 \n",
       "L 279.604253 139.161352 \n",
       "L 280.573161 138.095933 \n",
       "L 281.495804 136.773111 \n",
       "L 282.360855 135.121187 \n",
       "L 283.173324 133.149164 \n",
       "L 285.519669 126.378963 \n",
       "L 286.456873 124.857436 \n",
       "L 287.569154 124.299593 \n",
       "L 288.851685 124.733408 \n",
       "L 290.256214 125.929318 \n",
       "L 293.153545 129.066883 \n",
       "L 294.548066 130.415126 \n",
       "L 295.878264 131.420444 \n",
       "L 297.150394 132.101747 \n",
       "L 298.374971 132.50857 \n",
       "L 299.560688 132.68238 \n",
       "L 300.720225 132.696804 \n",
       "L 302.986632 132.393234 \n",
       "L 305.217167 131.858146 \n",
       "L 309.635184 130.497137 \n",
       "L 319.53162 127.141777 \n",
       "L 328.312593 124.023346 \n",
       "L 330.496047 123.034145 \n",
       "L 331.57899 122.356 \n",
       "L 332.652839 121.46101 \n",
       "L 333.708475 120.113074 \n",
       "L 334.737432 118.02584 \n",
       "L 335.725867 114.697226 \n",
       "L 336.660871 109.49029 \n",
       "L 337.547764 102.071648 \n",
       "L 338.431832 93.238407 \n",
       "L 339.416239 86.128005 \n",
       "L 340.589051 84.750887 \n",
       "L 341.914799 89.811814 \n",
       "L 344.554439 105.04092 \n",
       "L 345.76739 110.129104 \n",
       "L 346.925332 113.17268 \n",
       "L 348.051653 114.812108 \n",
       "L 349.16159 115.587256 \n",
       "L 350.263529 115.833194 \n",
       "L 351.361942 115.812615 \n",
       "L 353.55511 115.346054 \n",
       "L 356.842933 114.326695 \n",
       "L 366.703214 110.968921 \n",
       "L 378.746995 106.828219 \n",
       "L 394.060272 101.580823 \n",
       "L 409.35685 96.342647 \n",
       "L 419.184821 92.960045 \n",
       "L 424.641382 91.086478 \n",
       "L 424.641382 91.086478 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 136.060306 181.883287 \n",
       "L 163.831518 172.387244 \n",
       "L 168.256024 170.845525 \n",
       "L 172.703062 169.345758 \n",
       "L 177.127075 167.807559 \n",
       "L 183.786064 165.548119 \n",
       "L 189.321626 163.644619 \n",
       "L 203.710488 158.709976 \n",
       "L 240.173664 146.199791 \n",
       "L 246.799888 143.947655 \n",
       "L 257.824414 140.154822 \n",
       "L 264.430818 137.86081 \n",
       "L 269.884811 135.732325 \n",
       "L 272.012602 134.63217 \n",
       "L 273.052139 133.96388 \n",
       "L 274.063203 133.155356 \n",
       "L 275.046043 132.201219 \n",
       "L 275.977056 130.98196 \n",
       "L 276.854536 129.476579 \n",
       "L 277.674266 127.650395 \n",
       "L 279.165041 123.096543 \n",
       "L 279.91458 120.745044 \n",
       "L 280.764887 118.857038 \n",
       "L 281.781757 117.809914 \n",
       "L 282.985773 117.761122 \n",
       "L 284.361266 118.677916 \n",
       "L 287.306594 121.780096 \n",
       "L 288.756231 123.307425 \n",
       "L 290.134131 124.479774 \n",
       "L 291.453989 125.359542 \n",
       "L 292.711214 125.900352 \n",
       "L 293.921245 126.179402 \n",
       "L 295.10196 126.29621 \n",
       "L 297.388194 126.086943 \n",
       "L 299.625517 125.583395 \n",
       "L 301.84327 124.961029 \n",
       "L 308.449346 122.790782 \n",
       "L 312.846144 121.30822 \n",
       "L 319.433383 119.033117 \n",
       "L 324.919831 117.121312 \n",
       "L 327.11143 116.315129 \n",
       "L 329.292866 115.32883 \n",
       "L 330.374276 114.650269 \n",
       "L 331.445657 113.746635 \n",
       "L 332.500119 112.437224 \n",
       "L 333.523132 110.329877 \n",
       "L 334.505375 107.045567 \n",
       "L 335.427523 101.846602 \n",
       "L 336.29334 94.380839 \n",
       "L 337.149642 85.428131 \n",
       "L 338.102966 77.959648 \n",
       "L 339.261572 76.231861 \n",
       "L 340.597556 81.131396 \n",
       "L 343.291058 96.743802 \n",
       "L 344.523192 102.015818 \n",
       "L 345.693876 105.248281 \n",
       "L 346.826902 106.978729 \n",
       "L 347.940181 107.796105 \n",
       "L 349.04379 108.104222 \n",
       "L 350.142571 108.077177 \n",
       "L 352.335122 107.650448 \n",
       "L 357.810939 105.899091 \n",
       "L 364.380366 103.630583 \n",
       "L 372.039341 101.030074 \n",
       "L 376.415807 99.512304 \n",
       "L 382.975081 97.28265 \n",
       "L 390.627167 94.642534 \n",
       "L 397.179263 92.417164 \n",
       "L 402.639297 90.542849 \n",
       "L 412.464989 87.154646 \n",
       "L 416.824756 85.682921 \n",
       "L 416.824756 85.682921 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 128.345246 176.379124 \n",
       "L 136.103308 173.690802 \n",
       "L 140.561523 172.194615 \n",
       "L 144.990275 170.656683 \n",
       "L 149.437233 169.149247 \n",
       "L 153.871175 167.623871 \n",
       "L 163.853613 164.211068 \n",
       "L 170.490882 161.912101 \n",
       "L 174.938688 160.424842 \n",
       "L 179.354326 158.88052 \n",
       "L 183.793653 157.384533 \n",
       "L 188.214804 155.855784 \n",
       "L 194.862209 153.600892 \n",
       "L 199.273908 152.058653 \n",
       "L 203.708453 150.57048 \n",
       "L 208.119867 149.032695 \n",
       "L 212.547896 147.537831 \n",
       "L 216.962681 146.014276 \n",
       "L 223.595332 143.763272 \n",
       "L 228.000361 142.221774 \n",
       "L 232.424379 140.737872 \n",
       "L 236.830297 139.205346 \n",
       "L 241.245989 137.706814 \n",
       "L 247.854532 135.418989 \n",
       "L 252.271557 133.939079 \n",
       "L 256.668328 132.391019 \n",
       "L 259.967289 131.236354 \n",
       "L 263.247446 130.00613 \n",
       "L 266.466727 128.51132 \n",
       "L 268.533639 127.151385 \n",
       "L 269.521291 126.254834 \n",
       "L 270.464565 125.142849 \n",
       "L 271.362717 123.802382 \n",
       "L 272.183132 122.071265 \n",
       "L 272.949718 120.044562 \n",
       "L 275.13818 113.098585 \n",
       "L 276.044663 111.549675 \n",
       "L 277.158836 111.025899 \n",
       "L 278.472427 111.535787 \n",
       "L 279.919429 112.782697 \n",
       "L 282.926737 116.070504 \n",
       "L 284.366569 117.467678 \n",
       "L 285.735614 118.527919 \n",
       "L 287.030085 119.210588 \n",
       "L 288.277789 119.658086 \n",
       "L 289.47546 119.837331 \n",
       "L 291.793328 119.780145 \n",
       "L 295.155639 119.049094 \n",
       "L 302.873318 116.602156 \n",
       "L 310.556756 113.95039 \n",
       "L 314.949826 112.476742 \n",
       "L 318.23942 111.326981 \n",
       "L 325.912889 108.630182 \n",
       "L 328.090731 107.621208 \n",
       "L 329.172165 106.968258 \n",
       "L 330.240229 106.040759 \n",
       "L 331.293327 104.760986 \n",
       "L 332.312372 102.671811 \n",
       "L 333.286435 99.378763 \n",
       "L 334.198482 94.237625 \n",
       "L 335.042942 86.705759 \n",
       "L 335.871878 77.635 \n",
       "L 336.794429 69.847859 \n",
       "L 337.935414 67.717306 \n",
       "L 339.280752 72.475787 \n",
       "L 342.027134 88.425779 \n",
       "L 343.279952 93.93635 \n",
       "L 344.462944 97.298157 \n",
       "L 345.603791 99.173252 \n",
       "L 346.720243 100.013476 \n",
       "L 347.82572 100.371887 \n",
       "L 348.924999 100.371113 \n",
       "L 351.117044 99.961694 \n",
       "L 358.777952 97.449174 \n",
       "L 365.340376 95.227748 \n",
       "L 369.715805 93.693064 \n",
       "L 372.994599 92.586801 \n",
       "L 381.7359 89.602717 \n",
       "L 409.025378 80.22852 \n",
       "L 409.025378 80.22852 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 120.607111 170.836966 \n",
       "L 126.166953 168.946151 \n",
       "L 187.110764 148.080642 \n",
       "L 191.547632 146.596687 \n",
       "L 195.957081 145.057188 \n",
       "L 200.385903 143.56377 \n",
       "L 204.800202 142.040751 \n",
       "L 211.433859 139.792242 \n",
       "L 215.838567 138.253166 \n",
       "L 220.262674 136.768167 \n",
       "L 224.66884 135.239024 \n",
       "L 229.085283 133.742519 \n",
       "L 235.693918 131.458554 \n",
       "L 240.112078 129.980388 \n",
       "L 244.510542 128.443158 \n",
       "L 248.921889 126.953896 \n",
       "L 253.315471 125.405607 \n",
       "L 256.605592 124.225645 \n",
       "L 259.859576 122.899793 \n",
       "L 261.973122 121.778878 \n",
       "L 263.004924 121.107104 \n",
       "L 263.998509 120.265656 \n",
       "L 264.965717 119.298355 \n",
       "L 265.867923 118.031649 \n",
       "L 266.710968 116.478618 \n",
       "L 267.487686 114.594864 \n",
       "L 270.375334 105.564848 \n",
       "L 271.383267 104.545472 \n",
       "L 272.600105 104.531249 \n",
       "L 274.01779 105.529638 \n",
       "L 277.068932 108.772399 \n",
       "L 278.569364 110.362182 \n",
       "L 279.986475 111.575519 \n",
       "L 281.336885 112.482893 \n",
       "L 282.618523 113.056989 \n",
       "L 283.84176 113.338227 \n",
       "L 285.036461 113.480901 \n",
       "L 287.33298 113.276919 \n",
       "L 289.574081 112.777325 \n",
       "L 291.794383 112.167383 \n",
       "L 302.790626 108.531709 \n",
       "L 306.079951 107.38593 \n",
       "L 324.714024 100.924238 \n",
       "L 326.889778 99.91484 \n",
       "L 327.971559 99.290343 \n",
       "L 329.037854 98.366959 \n",
       "L 330.087558 97.076048 \n",
       "L 331.10442 95.034417 \n",
       "L 332.070629 91.7349 \n",
       "L 332.971686 86.616118 \n",
       "L 333.797894 79.080236 \n",
       "L 334.597566 69.83698 \n",
       "L 335.490678 61.787486 \n",
       "L 336.611335 59.230912 \n",
       "L 337.963964 63.824652 \n",
       "L 340.762625 80.091539 \n",
       "L 342.03721 85.858701 \n",
       "L 343.233116 89.365826 \n",
       "L 344.381677 91.347858 \n",
       "L 345.502197 92.270111 \n",
       "L 346.60906 92.624853 \n",
       "L 347.709314 92.685544 \n",
       "L 349.900797 92.271079 \n",
       "L 354.276366 90.872288 \n",
       "L 358.649401 89.426379 \n",
       "L 363.022814 87.888376 \n",
       "L 367.392436 86.434804 \n",
       "L 370.670457 85.298877 \n",
       "L 375.039602 83.797975 \n",
       "L 384.864467 80.440502 \n",
       "L 391.413566 78.173831 \n",
       "L 394.684111 77.071385 \n",
       "L 401.223051 74.860633 \n",
       "L 401.223051 74.860633 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "    <path clip-path=\"url(#p5ae19fdc2b)\" d=\"M 112.914117 165.351673 \n",
       "L 120.669603 162.672374 \n",
       "L 125.121775 161.173443 \n",
       "L 131.765731 158.879693 \n",
       "L 137.312455 156.987143 \n",
       "L 147.298648 153.594792 \n",
       "L 155.039371 150.921903 \n",
       "L 159.478473 149.425468 \n",
       "L 166.108504 147.135163 \n",
       "L 171.644841 145.252007 \n",
       "L 225.753472 126.76473 \n",
       "L 232.360269 124.489981 \n",
       "L 238.97425 122.247444 \n",
       "L 253.237391 117.183448 \n",
       "L 256.444085 115.680327 \n",
       "L 258.488914 114.28732 \n",
       "L 259.462849 113.379876 \n",
       "L 260.382114 112.231485 \n",
       "L 261.249089 110.842537 \n",
       "L 262.032017 109.06774 \n",
       "L 262.746739 106.959895 \n",
       "L 264.777002 99.830236 \n",
       "L 265.661657 98.291748 \n",
       "L 266.775644 97.785344 \n",
       "L 268.121216 98.367322 \n",
       "L 269.616176 99.693201 \n",
       "L 272.73165 103.122696 \n",
       "L 274.21371 104.557026 \n",
       "L 275.618806 105.657575 \n",
       "L 276.940492 106.374918 \n",
       "L 278.204466 106.825228 \n",
       "L 279.418244 107.035356 \n",
       "L 281.748962 106.982489 \n",
       "L 285.11905 106.268166 \n",
       "L 291.730805 104.16675 \n",
       "L 296.12548 102.706302 \n",
       "L 302.704069 100.442295 \n",
       "L 319.143405 94.828593 \n",
       "L 324.607779 92.790028 \n",
       "L 326.771793 91.606526 \n",
       "L 327.837899 90.712441 \n",
       "L 328.883957 89.405643 \n",
       "L 329.897989 87.392891 \n",
       "L 330.858187 84.118158 \n",
       "L 331.747635 78.99799 \n",
       "L 332.556944 71.473035 \n",
       "L 333.32842 62.08145 \n",
       "L 334.190315 53.737659 \n",
       "L 335.290496 50.801337 \n",
       "L 336.646681 55.16134 \n",
       "L 339.497759 71.753414 \n",
       "L 340.794731 77.773965 \n",
       "L 342.004482 81.450952 \n",
       "L 343.160659 83.519924 \n",
       "L 344.285729 84.532787 \n",
       "L 345.394191 84.905747 \n",
       "L 346.495199 84.985845 \n",
       "L 348.686242 84.575917 \n",
       "L 350.873198 83.922444 \n",
       "L 354.151882 82.84375 \n",
       "L 358.521828 81.37605 \n",
       "L 362.890987 79.872213 \n",
       "L 371.624665 76.894548 \n",
       "L 393.439533 69.412513 \n",
       "L 393.439533 69.412513 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-opacity:0.5;\"/>\n",
       "   </g>\n",
       "  </g>\n",
       " </g>\n",
       " <defs>\n",
       "  <clipPath id=\"p5ae19fdc2b\">\n",
       "   <rect height=\"334.8\" width=\"669.6\" x=\"7.2\" y=\"7.2\"/>\n",
       "  </clipPath>\n",
       " </defs>\n",
       "</svg>\n"
      ],
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10ce7a2d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Set up grid and differentiation matrix:\n",
    "N = 256; dt = 0.4/N**2; x = (2*pi/N)*arange(-N/2,N/2);\n",
    "A, B = 25.0, 16.0\n",
    "u = 3*A**2/cosh(0.5*A*(x+2))**2 + 3*B**2/cosh(0.5*B*(x+1))**2\n",
    "v = fft(u); \n",
    "k = zeros(N); k[0:N/2] = arange(0,N/2); k[N/2+1:] = arange(-N/2+1,0,1)\n",
    "ik3 = 1j*k**3\n",
    "\n",
    "# Time-stepping by Runge-Kutta\n",
    "tmax = 0.006; nplt = int(round((tmax/25)/dt))\n",
    "nmax = int(round(tmax/dt))\n",
    "udata = []; udata.append(list(zip(x, u)))\n",
    "tdata = [0.0]\n",
    "for n in range(1,nmax+1):\n",
    "    t = n*dt; g = -0.5j*dt*k\n",
    "    E = exp(dt*ik3/2); E2 = E**2\n",
    "    a = g * fft(real(ifft( v         ))**2)\n",
    "    b = g * fft(real(ifft( E*(v+a/2) ))**2)\n",
    "    c = g * fft(real(ifft( E*v+b/2   ))**2)\n",
    "    d = g * fft(real(ifft( E2*v+E*c  ))**2)\n",
    "    v = E2*v + (E2*a + 2*E*(b+c) + d)/6\n",
    "    if n%nplt == 0:\n",
    "        u = real(ifft(v))\n",
    "        udata.append(list(zip(x, u)))\n",
    "        tdata.append(t);\n",
    "\n",
    "fig = figure(figsize=(12,6))\n",
    "ax = fig.gca(projection='3d')\n",
    "poly = LineCollection(udata)\n",
    "poly.set_alpha(0.5)\n",
    "ax.add_collection3d(poly, zs=tdata, zdir='y')\n",
    "ax.set_xlabel('x')\n",
    "ax.set_xlim3d(-pi, pi)\n",
    "ax.set_ylabel('t')\n",
    "ax.set_ylim3d(0, tmax)\n",
    "ax.set_zlabel('u')\n",
    "ax.set_zlim3d(0, 2000)\n",
    "ax.view_init(80,-125);"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [default]",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}
